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The average current and the shot noise at correlated se- 
quential tunneling via two localized sites are studied. At zero 
temperature the Fano factor averaged over the positions and 
energies of sites is shown to be 0.707. The noise dependence 
on temperature and frequency is analyzed numerically. 
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Shot noise in mesoscopic structures has been the sub- 
ject of thorough studies in the recent past.EFB In partic- 
ular, a theory of shot noise at tunneling of singlcuclcc- 
trons correlated due^ta-. Coulomb blockade effcctsa has 
been well develop ec&119 and verified experimentallyxil 
Recently the first attempts have been jaadeOtj to ex- 
tend this theory to hopping transportlU which can be 
formally considered as a special case of correlated single- 
electron tunneling. 

In a typical hopping situation there is a considerable 
1/f contribution at low frequencies (see Ref.li3 and ref- 
erences therein), so one can discuss the shot noise only 
at sufficiently high frequencies. The 1/f noise at hop- 
ping is mainly due to electron-electron interaction: the 
slowly evolving trapped charge configurations can signif- 
icantly affect the current through nearby channels. The 
1// component is absent-at hopping through noninter- 
acting ID chains of siteslia (while the slow fluctuations 
of the chain parameters due to external traps can restore 
this component). In the present paper we consider hop- 
ping through very short chains which are just pairs of 
sites, and assume that the parameters of these pairs do 
not fluctuate in time, so that the noise does not have 1/f 
contribution. 

For two-site hopping we basically follow the model in- 
troduced by Glazman and Matveev.EJ The only diffcr- 



correlation between 
iin this respect our 



ence is that we take into account t 
tunneling events neglected in Ref.t 
model is closer to the model of Ref.lHI) . Using the meth- 
ods developed in Ref.H we calculate the current / and 
the current spectral density Si(u) for an individual two- 
site channel. The summation over many parallel cliaa- 
nels with random parameters is done similar to Ref.EJl 
The main object of our study is the Fano factor F (the 
low frequency noise normalized by the Schottky value 
Sj = 2el). We will show that at zero temperature the 
Fano factors for individual two-site channels range from 
5/14 to 1, while after averaging we get F — 0.707 (the 
similarr-Broblem for one-site channels has been considered 
in Refill with the result F = 3/4). For a finite temper- 
ature T the Fano factor can be calculated numerically; 
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FIG. 1. Schematic of two-site tunneling channel. 

after the averaging we obtain F as a function of the ratio 
T/eV where V is the voltage between electrodes. 

The schematic of a two-site channel is shown in Fig. 
[]]. The thickness d of an insulating layer between two 
metallic electrodes is assumed to be much greater than 
the electron localization radius a, and we use the model 
of sequential (incoherent) hops of single electrons. The 
lengths of the left and right hops are x\ and £3, respec- 
tively, while the hop between two sites has the length 
r 2 = ( x 2 + U 2 ) 1 ^ 2 where X2 — d — X\ — X3 and y is the 
shift of site positions in the plane parallel to electrodes. 
Each site can be occupied by at most one electron and 
the effect of electron spin is neglected (the case of double- 
degeneracy due to spin will be discussed later). The tun- 
neling rates from electrodes to empty sites (tunneling to 
nearest neighbor only, see Fig. pi) are assumed to be 



T+ = T cxp(-2xi/a)/(~ey + £1), 

= r cx P (-2 a;3 / a )/(£ 2 ), 



(i) 



where superscripts indicate the direction of tunneling, e\ 
and £2 are the site energies counted from the Fermi level 
of the right electrode, and f(e) — [1 + exp(£/T)] _1 is the 
Fermi function. Similarly, the rates of tunneling from 
occupied sites to neighboring electrodes are 



T- = r exp(-2xi/a)/(ey - £ X ), 
r+ = r exp(-2z 3 /a)/(-£2). 



(2) 



Notice that we have neglected the Coulomb interaction 
of electrons on different sites (energies £1,2 do not depend 
on the occupation of neighboring site). The rate of in- 
elastic tunneling between the sites depends on the energy 
difference Ae — e\ — £2, and for |Ae| much smaller than 
(where s is the sound velocity) can be calculated 



T2 = cvFq exp(— 2r2/a) 



±Ae 



1 - exp(=FA£/T) ' 



(3) 
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FIG. 2. Graphical representation of the master equation 
for zero temperature. 



where the dimensional factor a describes the relative 
strength of phonon-assisted tunneling compared to "res- 
onant" tunneling assumed in Eqs. ([!])- (||). 

Let us start with zero temperature case. Then the 
transport is possible only if eV > e\ > £2 > 0, and elec- 
trons move only in one direction, = T 2 = r 3 = 
(for simplicity we omit the superscript "+", I\ = 
The kinetic ( "master" ) equation in this case can be rep- 
resented graphically by Fig. |^. The configuration space 
consists of four charge states of the two-site system which 
are denoted as 00, 01, 10, and 11, while arrows show 
transitions between them. The graphical representation 
of the master equation in a relatively small configura- 
tion space is a very convenient tool and often allows 
straightforward calculation of the average current and 
zero-frequency spectral density (see, e.g., Refill). 

The basic idea of the methods is to consider the ran- 
dom "travel" of the system state within the configura- 
tion space and divide the duration of this stochastic pro- 
cess into blocks which start and end in a specific charge 
state. Because of the Markovian property of the process 
(absence of memory) these blocks are mutually uncorre- 
cted, so the averaging over the blocks is rather simple. 
In the case of Fig. |] let us choose the charge state 01 as 
the block divider. Then there are two types of blocks: 
01 -» 00 -> 10 -> 01 (type 1) and 01 -> 11 ->• 10 -> 01 
(type 2), while the blocks are additionally characterized 
by the time spent in each charge state. n 

The average current can be calculatedQ as 



/ = efe/r, 



(4) 



where f is the average block duration and k is the av- 
erage number of electrons transferred between electrodes 
per block (the averaging is taken over a large number 
of blocks). To calculate these average magnitudes let us 
notice that the blocks of type 1 and type 2 have proba- 
bilities 



pi = r 3 /(r x + r 3 ), P2 = iy^ + r 3 



(5) 



and the average durations t\ and t 2 of the blocks of each 
type can be calculated as 



n = (r 1 + r 3 )- 1 + rr 1 + r 2 - 1 , 
T2 = (r 1 + r3)- 1 + r 3 - 1 + r 2 - 1 



(6) 



(notice that the average waiting time (Fi + r^ -1 of the 
first hop is equal for both types). Taking into account 
that each block corresponds to the transfer of one elec- 
tron, k\ = k<i = k = 1, and calculating the average block 
duration 

T=p 1 T\+p 2 T2, (7) 

we finally obtain the formula for the average current @ 



1 1 + ri/r 3 + r 3 /r x 
r 2 ri + r 3 



(8) 



This equation differs from the result of RefB because 
the correlation between the occupations of two sites was 
neglected in RefJiS. The correct equation for the current 
which coincides with Eq. (||) was obtained later in Ref.tZI. 

The same method as above can be used for the cal- 
culation of the low-frequency limit 5/(0) of the curre 
spectral density which is given by the general equatio 
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(0) = (2/r) (e 2 k 2 + I 2 t 2 - 2elkrj 



(9) 



(averaging is again over blocks) which in our case at zero 
temperature simplifies to 

5/(0) = 2el [P/t 2 ) - lj . (10) 
So, besides Eqs. (§)-(0) we also need to calculate t 2 : 

^ = Pi^ + P2^f, (11) 

where because of Poissonian statistics of each tunneling 
event we have 



2 , p-2 

/■j - 1 1 ! -,- j. 3) f r 3 2 + r 2 2 . 



T 2 -w 2 = (ri + r 3 )- 2 + r^ 
?-T2 2 = (r 1 + r r ^ 2 



(12) 



Combining these equations we finally obtain 



1 (l + R + R- 1 ) 2 -4 



5/(0) - 2el ( ^2 



(ri + r 3 )2 



1 l + R + R' 1 

ri + r 3 



(13) 



where R = Ti/r 3 . Analyzing the Fano factor F = 
5/(0)/2eJ one can see that the uniform case, ri = T2 = 
r 3 , provides F = 9/25 which is not the minimum pos- 
sible value. The minimum Fano factor is achieved at 
Ti = T 3 = 5/6 x T 2 and equal to F mm = 5/14 (it is still 
noticeably larger than the naive estimate F = 1/3). The 
maximal value F max = 1 is obviously achieved when one 
of the rates Ti isanuch smaller than two other rates. 

Following Ref.li3 let us assume many two-site channels 
"in parallel" and find the total current 1-% and spectral 
density 5/^(0) integrating over channels with different 
site positions x\, X3, y and different energies E\ and £ 2 . 
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Assuming a sufficiently thick insulating layer we may ap- 
proximate the distance between sites as r-i — x2 + y 2 /2x2, 
where X2 corresponds to the channel with the maximum 
current. For such a channel y — and ri — T2 = T3 = 
T exp(-2d/3a)(aeV)^ 3 [see Eq. (§) and Eqs. 
which gives X2 = (d/3) + (a/3) \n(aeV). At zero temper- 
ature the total current can be calculated as 



Iv=n z A / da / d6 / lirydy 

J — 00 J — 00 Jo 

r-eV 

x dAe (eV- Ae) I(^,^,y, As), (14) 
Jo 

where n is the density of states, A is the area [A ^> d 2 , 
A 3> n~ 2 a~ 3 d~ 1 (eV)~ 2 ], the x-positions of the sites are 
measured from the optimal values, £j — Xi — Xi, x% = 
X3 = (d — X2)/2 (the integration is extended to infinity 
since d 3> a), and the current / is given by Eq. (||). 

Using the relation I(£,i,£,3,y, Ae) = exp(— 25 /a) I(£i — 
S, £3 — 6, 0, Ae) where 5 — y 2 /6x2, it is easy to show that 
the integration L 2irydy gives the factor 3irax2- Cal- 
culating the integral over Ae analytically and integrals 
over £1 and £3 numerically, we get the result 

Is = 5.237 eT n 2 Aa 3 X2 exp(-2d/3a) a 1/3 (eV) 7/3 . (15) 



Notice that the scaling Is oc y 7 / 3 _is the same as in the 
model which neglects correlations £3 

A similar sum over different two-site channels can be 
calculated for the current spectral density at zero fre- 
quency [just replacing / in Eq. ( fhij ) with 57(0) given by 
Eq. (]13|)1. Integrating analytically over y and Ae and 
numerically over the two remaining variables we obtain 
the following average Fano factor at zero temperature: 



F = 5/, s (0)/2e/ s = 0.7074. 



(16) 



Let us now consider the finite-temperature case. Our 
method for calculation of I and 5/(0) based on the anal- 
ysis of blocks can still be easily applied if T^" ^ while 
IY = IY = (this situation occurs when s\ and £2 are 
well inside the energy strip defined by the Fermi levels of 
the electrodes). In this case the current I and the specu 
tral density 57(0) are given by Eqs. (Q) and (||), whereEj 

t = pin + P2T2 + P3T3, k = k 2 =p!+p 2 , 

pi = r 3"/ r s, P2 = rf/r E , p 3 = r^/rs, 

n = i/r E + i/r+ + i/r+ n = i/r E + i/r+, 

t- 2 = i/r E + i/r+ + i/r+ r s = r+ + r+ + r 2 - , 

t 2 = pirl + p 2 rl + p 3 r 3 2 , kr = pin + P2T2, 

?-n 2 = (r s )- 2 + (r+r 2 + (r+r 2 , 



^-T2 2 = (r s )- 2 + (r 3 +)- 2 + (r+)- 2 . 
r|-T3 2 = (r s )- 2 + (r+)- 2 . 



(17) 



In the case when all 1^ are non-zero, it is more nat- 
ural to use the general master-equation formalism for 
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FIG. 3. The Fano factor F averaged over two-site (thick 
solid line) and one-site (thin solid line) channels, as a function 
of temperature. Dotted lines show F(T) for direct tunnel- 
ing and Ohmic conduction. Dashed line shows the averaged 
two-site current I(T) normalized by 1(0). 



the average currenta and spectral density.^ We have de- 
veloped a numerical code and integrated over different 
two-site channels in the same way as above, just us- 
ing the numerical results for / and 5/(0) instead of 
Eqs. (g) and @. (One more difference from the zero- 
temperature case is the separate integration over ei and 
£2 ) The dashed line in Fig. |3| shows the dependence 
of the ratio g = Is(T) / /s(0) on the normalized tem- 
perature T/eV. [We neglect the weak temperature de- 
pendence of X2 and actually calculate the dependence 
of the numerical factor in Eq. (|l5|).] The asymptote at 
T > eV is g = 21.7 (T/eF) 4 / 3 , so the conductance G is 
equal to G = 113.6 e 2 r n 2 Aa 3 x 2 exp(-2d/3a) a^ 3 T 4 / 3 , 
where X2 can be approximated as (d/3) + (a/3)ln(aT) 
(the scalingJT 4 / 3 is the same as in the model without 
correlationsES) . 

The Fano factor averaged over different channels, as a 
function of T/eV is shown in Fig. || by the thick solid line. 
The low-temperature value is given by Eq. (|l6|), while the 
high-temperature asymptote, F = 2T/eV (lower dotted 
line) can be easily derived from the Nyquist formula. 

It is interesting to compare the temperature depen- 
dence of F for two-site and one-site channels. In the 
latter case we still use Eqs. (|]) and (§) for the tunnel- 
ing rates similar to Ref.lla. The thin solid line in Fig. 
H shows the average Fano factor for one-site channels, 
as a function of T/eV (this.-c.urve in other coordinates 
has been calculated in Ref.E£l). The low temperature 
value is F = 3/4, while the high-temperature asymp- 
tote, F = 2T/eV, is the same as for two-site channels 
and direct-tunneling case. (The result for the direct 
tunneling, F = coth(el//2T), is shown for comparison 
by the upper dotted line.) Obviously, with the increase 
of the number N of sites in the channel the average Fano 
factor decreases. However, its dependence on N seems to 
saturate rapidly, as indicated by the small difference be- 
tween the results for one-site and two-site channels. So, 
even for large TV one should expect the dependence F(T) 
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FIG. 4. Solid line: the frequency dependence of the normal- 
ized current spectral density Si (u))/2el averaged over two-site 
channels. Dashed lines correspond to particular channels (see 
text). The dotted line shows Si(ui)/2el averaged over one-site 
channels. 



to deviate significantly at T < eV from the result for an 
Ohmic conductor, F = 2T/eV. This can be explained 
by the fact that the ID chains of sites with "soft" (not 
strong) bottlenecks still give considerable contribution to 
the total current, while the Fano factor for such-chains is 
comparable to unity. (The situation is differentt2l for 2D 
or 3D hopping because the percolation cluster does not 
have bottlenecks at the size scale much larger than the 
correlation length of the cluster. As a consequence, for 
sufficiently large samples we expect F — 2T/eV.) 

So far we have discussed only the current spectral 
density at zero frequency. Our computer code can also 
treat the finite-frequency case. At finite frequency w it 
is necessary to specify where the current is measured. 
We have considered the current in the electrodes and 
assumed natural electrostatics when the electron hop 
through ith gap transfers charge qi = exi/d in the elec- 
trodes. The spectral density at finite u) depends not only 
on tunneling rates Tf but also on and thus on the 
positions of sites. For averaging over the two-site chan- 
nels we have used the approximation qi ~ exi/d ~ e/3. 
The solid line in Fig. ^ shows the frequency dependence 
5/ ! s(^)/2e/s for T = 0. The frequency scale is de- 
termined by tunneling rates, so for normalization we 
have used r„ = Fo cxp(— 2d/3a)(aeV) 1 ^ 3 . For compar- 
ison, we also show Si(w)/2el for the uniform channel, 
ri = T2 = T3 (lower dashed line) and for two nonuni- 
form channels: T\ = T2 = O.IT3 (middle dashed line) 
and lOTi = T2 = O.IT3 (upper dashed line) at Ae = eV. 
Notice that the solid line (averaged noise) has a finite 
slope at iij = (even though the slope is zero for each 
individual channel) and approaches the high-w asymp- 
tote of 1/3 as lo^ 1 / 2 . The dotted line in Fig. |J shows 
the ratio Si(ui)/2el averaged over one-site channels [then 
r n = Fo exp(— d/a)], which can be calculated analyti- 
cally: 1/2 + [{uj 2 T- 2 + A) 1 ' 2 - 2] T 2 Jlj 2 . 

Finally, let us briefly consider the effect of electron spin 
using a simple model. Assuming the double degeneracy 



due to spin (but still allowing at most one electron per 
site), we should double the tunneling rates r+ and T^" 
[see Eq. (p|)], while leaving all other rates unchanged. At 
zero temperature this will lead to a trivial extra factor 
Eq. @ and a very small change of X2 , while the 
average Fano factor given by Eq. ( |l6|) does not change. 
The calculations at finite temperature show that F is a 
little larger in the case of double degeneracy compared 
to the spinless case, however, the difference is so small 
that the corresponding curves in Fig. ^cannot be dis- 
tinguished. The maximum difference AF « 5 x 10~ 4 is 
achieved at T/eV « 0.3, while at T > eV the differ- 
ence approaches zero. A similar very weak dependence 
on the spin degeneracy for one-site channels has been re- 
ported in Ref.ta (we have found the maximum difference 
AF = 1.12 x 1CT 3 at T/eV « 0.33). 

In conclusion, we have studied the shot noise of two- 
site hopping channels. The different average Fano factor 
and different frequency dependence of the noise in com- 
parison with one-site channels and direct tunneling can in 
principle be verified experimentally (using the difference 
of the temperature dependence of the average current). 
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